Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders.

Purpose: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). Methods: We coupled phenotyping with exome or genome sequencing of 467 pedigrees with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. Results: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. Conclusion: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

The recurrent deep intronic pseudoexon-inducing variant COL6A1 c.930+189C>T results in a consistently severe phenotype of COL6-related dystrophy: Towards clinical trial readiness for splice-modulating therapy.

Collagen VI-related dystrophies (COL6-RDs) manifest with a spectrum of clinical phenotypes, ranging from Ullrich congenital muscular dystrophy (UCMD), presenting with prominent congenital symptoms and characterised by progressive muscle weakness, joint contractures and respiratory insufficiency, to Bethlem muscular dystrophy, with milder symptoms typically recognised later and at times resembling a limb girdle muscular dystrophy, and intermediate phenotypes falling between UCMD and Bethlem muscular dystrophy. Despite clinical and immunohistochemical features highly suggestive of COL6-RD, some patients had remained without an identified causative variant in COL6A1, COL6A2 or COL6A3. With combined muscle RNA-sequencing and whole-genome sequencing we uncovered a recurrent, de novo deep intronic variant in intron 11 of COL6A1 (c.930+189C>T) that leads to a dominantly acting in-frame pseudoexon insertion. We subsequently identified and have characterised an international cohort of forty-four patients with this COL6A1 intron 11 causative variant, one of the most common recurrent causative variants in the collagen VI genes. Patients manifest a consistently severe phenotype characterised by a paucity of early symptoms followed by an accelerated progression to a severe form of UCMD, except for one patient with somatic mosaicism for this COL6A1 intron 11 variant who manifests a milder phenotype consistent with Bethlem muscular dystrophy. Characterisation of this individual provides a robust validation for the development of our pseudoexon skipping therapy. We have previously shown that splice-modulating antisense oligomers applied in vitro effectively decreased the abundance of the mutant pseudoexon-containing COL6A1 transcripts to levels comparable to the in vivo scenario of the somatic mosaicism shown here, indicating that this therapeutic approach carries significant translational promise for ameliorating the severe form of UCMD caused by this common recurrent COL6A1 causative variant to a Bethlem muscular dystrophy phenotype.

Genome and RNA sequencing boost neuromuscular diagnoses to 62% from 34% with exome sequencing alone.

OBJECTIVE: Most families with heritable neuromuscular disorders do not receive a molecular diagnosis. Here we evaluate diagnostic utility of exome, genome, RNA sequencing, and protein studies and provide evidence-based recommendations for their integration into practice. METHODS: In total, 247 families with suspected monogenic neuromuscular disorders who remained without a genetic diagnosis after standard diagnostic investigations underwent research-led massively parallel sequencing: neuromuscular disorder gene panel, exome, genome, and/or RNA sequencing to identify causal variants. Protein and RNA studies were also deployed when required. RESULTS: Integration of exome sequencing and auxiliary genome, RNA and/or protein studies identified causal or likely causal variants in 62% (152 out of 247) of families. Exome sequencing alone informed 55% (83 out of 152) of diagnoses, with remaining diagnoses (45%; 69 out of 152) requiring genome sequencing, RNA and/or protein studies to identify variants and/or support pathogenicity. Arrestingly, novel disease genes accounted for <4% (6 out of 152) of diagnoses while 36.2% of solved families (55 out of 152) harbored at least one splice-altering or structural variant in a known neuromuscular disorder gene. We posit that contemporary neuromuscular disorder gene-panel sequencing could likely provide 66% (100 out of 152) of our diagnoses today. INTERPRETATION: Our results emphasize thorough clinical phenotyping to enable deep scrutiny of all rare genetic variation in phenotypically consistent genes. Post-exome auxiliary investigations extended our diagnostic yield by 81% overall (34-62%). We present a diagnostic algorithm that details deployment of genomic and auxiliary investigations to obtain these diagnoses today most effectively. We hope this provides a practical guide for clinicians as they gain greater access to clinical genome and transcriptome sequencing.

Flavones provide resistance to DUX4-induced toxicity via an mTor-independent mechanism.

Facioscapulohumeral muscular dystrophy (FSHD) is among the most common of the muscular dystrophies, affecting nearly 1 in 8000 individuals, and is a cause of profound disability. Genetically, FSHD is linked to the contraction and/or epigenetic de-repression of the D4Z4 repeat array on chromosome 4, thereby allowing expression of the DUX4 gene in skeletal muscle. If the DUX4 transcript incorporates a stabilizing polyadenylation site the myotoxic DUX4 protein will be synthesized, resulting in muscle wasting. The mechanism of toxicity remains unclear, as many DUX4-induced cytopathologies have been described, however cell death does primarily occur through caspase 3/7-dependent apoptosis. To date, most FSHD therapeutic development has focused on molecular methods targeting DUX4 expression or the DUX4 transcript, while therapies targeting processes downstream of DUX4 activity have received less attention. Several studies have demonstrated that inhibition of multiple signal transduction pathways can ameliorate DUX4-induced toxicity, and thus compounds targeting these pathways have the potential to be developed into FSHD therapeutics. To this end, we have screened a group of small molecules curated based on their reported activity in relevant pathways and/or structural relationships with known toxicity-modulating molecules. We have identified a panel of five compounds that function downstream of DUX4 activity to inhibit DUX4-induced toxicity. Unexpectedly, this effect was mediated through an mTor-independent mechanism that preserved expression of ULK1 and correlated with an increase in a marker of active cellular autophagy. This identifies these flavones as compounds of interest for therapeutic development, and potentially identifies the autophagy pathway as a target for therapeutics.

Deep Mutational Scanning in Disease-related Genes with Saturation Mutagenesis-Reinforced Functional Assays (SMuRF).

Interpretation of disease-causing genetic variants remains a challenge in human genetics. Current costs and complexity of deep mutational scanning methods hamper crowd-sourcing approaches toward genome-wide resolution of variants in disease-related genes. Our framework, Saturation Mutagenesis-Reinforced Functional assays (SMuRF), addresses these issues by offering simple and cost-effective saturation mutagenesis, as well as streamlining functional assays to enhance the interpretation of unresolved variants. Applying SMuRF to neuromuscular disease genes FKRP and LARGE1, we generated functional scores for over 99.8% of all possible coding single nucleotide variants and resolved 310 clinically reported variants of uncertain significance with high confidence, enhancing clinical variant interpretation in dystroglycanopathies. SMuRF also demonstrates utility in predicting disease severity, resolving critical structural regions, and providing training datasets for the development of computational predictors. Our approach opens new directions for enabling variant-to-function insights for disease genes in a manner that is broadly useful for crowd-sourcing implementation across standard research laboratories.

Profilin1 is required for prevention of mitotic catastrophe in murine and human glomerular diseases.

The progression of proteinuric kidney diseases is associated with podocyte loss, but the mechanisms underlying this process remain unclear. Podocytes reenter the cell cycle to repair double-stranded DNA breaks. However, unsuccessful repair can result in podocytes crossing the G1/S checkpoint and undergoing abortive cytokinesis. In this study, we identified Pfn1 as indispensable in maintaining glomerular integrity - its tissue-specific loss in mouse podocytes resulted in severe proteinuria and kidney failure. Our results suggest that this phenotype is due to podocyte mitotic catastrophe (MC), characterized histologically and ultrastructurally by abundant multinucleated cells, irregular nuclei, and mitotic spindles. Podocyte cell cycle reentry was identified using FUCCI2aR mice, and we observed altered expression of cell-cycle associated proteins, such as p21, p53, cyclin B1, and cyclin D1. Podocyte-specific translating ribosome affinity purification and RNA-Seq revealed the downregulation of ribosomal RNA-processing 8 (Rrp8). Overexpression of Rrp8 in Pfn1-KO podocytes partially rescued the phenotype in vitro. Clinical and ultrastructural tomographic analysis of patients with diverse proteinuric kidney diseases further validated the presence of MC podocytes and reduction in podocyte PFN1 expression within kidney tissues. These results suggest that profilin1 is essential in regulating the podocyte cell cycle and its disruption leads to MC and subsequent podocyte loss.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality.

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

Death after High-Dose rAAV9 Gene Therapy in a Patient with Duchenne's Muscular Dystrophy.

We treated a 27-year-old patient with Duchenne's muscular dystrophy (DMD) with recombinant adeno-associated virus (rAAV) serotype 9 containing dSaCas9 (i.e., "dead" Staphylococcus aureus Cas9, in which the Cas9 nuclease activity has been inactivated) fused to VP64; this transgene was designed to up-regulate cortical dystrophin as a custom CRISPR-transactivator therapy. The dose of rAAV used was 1×1014 vector genomes per kilogram of body weight. Mild cardiac dysfunction and pericardial effusion developed, followed by acute respiratory distress syndrome (ARDS) and cardiac arrest 6 days after transgene treatment; the patient died 2 days later. A postmortem examination showed severe diffuse alveolar damage. Expression of transgene in the liver was minimal, and there was no evidence of AAV serotype 9 antibodies or effector T-cell reactivity in the organs. These findings indicate that an innate immune reaction caused ARDS in a patient with advanced DMD treated with high-dose rAAV gene therapy. (Funded by Cure Rare Disease.).

Development of Assays to Measure GNE Gene Potency and Gene Replacement in Skeletal Muscle.

BACKGROUND: GNE myopathy (GNEM) is a severe muscle disease caused by mutations in the UDP-GlcNAc-2-epimerase/ManNAc-6-kinase (GNE) gene, which encodes a bifunctional enzyme required for sialic acid (Sia) biosynthesis. OBJECTIVE: To develop assays to demonstrate the potency of AAV gene therapy vectors in making Sia and to define the dose required for replacement of endogenous mouse Gne gene expression with human GNE in skeletal muscles. METHODS: A MyoD-inducible Gne-deficient cell line, Lec3MyoDI, and a GNE-deficient human muscle cell line, were made and tested to define the potency of various AAV vectors to increase binding of Sia-specific lectins, including MAA and SNA. qPCR and qRT-PCR methods were used to quantify AAV biodistribution and GNE gene expression after intravenous delivery of AAV vectors designed with different promoters in wild-type mice. RESULTS: Lec3 cells showed a strong deficit in MAA binding, while GNE-/-MB135 cells did not. Overexpressing GNE in Lec3 and Lec3MyoDI cells by AAV infection stimulated MAA binding in a dose-dependent manner. Use of a constitutive promoter, CMV, showed higher induction of MAA binding than use of muscle-specific promoters (MCK, MHCK7). rAAVrh74.CMV.GNE stimulated human GNE expression in muscles at levels equivalent to endogenous mouse Gne at a dose of 1×1013vg/kg, while AAVs with muscle-specific promoters required higher doses. AAV biodistribution in skeletal muscles trended higher when CMV was used as the promoter, and this correlated with increased sialylation of its viral capsid. CONCLUSIONS: Lec3 and Lec3MyoDI cells work well to assay the potency of AAV vectors in making Sia. Systemic delivery of rAAVrh74.CMV.GNE can deliver GNE gene replacement to skeletal muscles at doses that do not overwhelm non-muscle tissues, suggesting that AAV vectors that drive constitutive organ expression could be used to treat GNEM.

Noncoding variants alter GATA2 expression in rhombomere 4 motor neurons and cause dominant hereditary congenital facial paresis.

Hereditary congenital facial paresis type 1 (HCFP1) is an autosomal dominant disorder of absent or limited facial movement that maps to chromosome 3q21-q22 and is hypothesized to result from facial branchial motor neuron (FBMN) maldevelopment. In the present study, we report that HCFP1 results from heterozygous duplications within a neuron-specific GATA2 regulatory region that includes two enhancers and one silencer, and from noncoding single-nucleotide variants (SNVs) within the silencer. Some SNVs impair binding of NR2F1 to the silencer in vitro and in vivo and attenuate in vivo enhancer reporter expression in FBMNs. Gata2 and its effector Gata3 are essential for inner-ear efferent neuron (IEE) but not FBMN development. A humanized HCFP1 mouse model extends Gata2 expression, favors the formation of IEEs over FBMNs and is rescued by conditional loss of Gata3. These findings highlight the importance of temporal gene regulation in development and of noncoding variation in rare mendelian disease.

Multi-omics identifies large mitoribosomal subunit instability caused by pathogenic MRPL39 variants as a cause of pediatric onset mitochondrial disease.

MRPL39 encodes one of 52 proteins comprising the large subunit of the mitochondrial ribosome (mitoribosome). In conjunction with 30 proteins in the small subunit, the mitoribosome synthesizes the 13 subunits of the mitochondrial oxidative phosphorylation (OXPHOS) system encoded by mitochondrial Deoxyribonucleic acid (DNA). We used multi-omics and gene matching to identify three unrelated individuals with biallelic variants in MRPL39 presenting with multisystem diseases with severity ranging from lethal, infantile-onset (Leigh syndrome spectrum) to milder with survival into adulthood. Clinical exome sequencing of known disease genes failed to diagnose these patients; however quantitative proteomics identified a specific decrease in the abundance of large but not small mitoribosomal subunits in fibroblasts from the two patients with severe phenotype. Re-analysis of exome sequencing led to the identification of candidate single heterozygous variants in mitoribosomal genes MRPL39 (both patients) and MRPL15. Genome sequencing identified a shared deep intronic MRPL39 variant predicted to generate a cryptic exon, with transcriptomics and targeted studies providing further functional evidence for causation. The patient with the milder disease was homozygous for a missense variant identified through trio exome sequencing. Our study highlights the utility of quantitative proteomics in detecting protein signatures and in characterizing gene-disease associations in exome-unsolved patients. We describe Relative Complex Abundance analysis of proteomics data, a sensitive method that can identify defects in OXPHOS disorders to a similar or greater sensitivity to the traditional enzymology. Relative Complex Abundance has potential utility for functional validation or prioritization in many hundreds of inherited rare diseases where protein complex assembly is disrupted.

Estimating the Prevalence of LAMA2 Congenital Muscular Dystrophy using Population Genetic Databases.

BACKGROUND: Recessive pathogenic variants in LAMA2 resulting in complete or partial loss of laminin α2 protein cause congenital muscular dystrophy (LAMA2 CMD). The prevalence of LAMA2 CMD has been estimated by epidemiological studies to lie between 1.36-20 cases per million. However, prevalence estimates from epidemiological studies are vulnerable to inaccuracies owing to challenges with studying rare diseases. Population genetic databases offer an alternative method for estimating prevalence. OBJECTIVE: We aim to use population allele frequency data for reported and predicted pathogenic variants to estimate the birth prevalence of LAMA2 CMD. METHODS: A list of reported pathogenic LAMA2 variants was compiled from public databases, and supplemented with predicted loss of function (LoF) variants in the Genome Aggregation Database (gnomAD). gnomAD allele frequencies for 273 reported pathogenic and predicted LoF LAMA2 variants were used to calculate disease prevalence using a Bayesian methodology. RESULTS: The world-wide birth prevalence of LAMA2 CMD was estimated to be 8.3 per million (95% confidence interval (CI) 6.27 -10.5 per million). The prevalence estimates for each population in gnomAD varied, ranging from 1.79 per million in East Asians (95% CI 0.63 -3.36) to 10.1 per million in Europeans (95% CI 6.74 -13.9). These estimates were generally consistent with those from epidemiological studies, where available. CONCLUSIONS: We provide robust world-wide and population-specific birth prevalence estimates for LAMA2 CMD, including for non-European populations in which LAMA2 CMD prevalence hadn't been studied. This work will inform the design and prioritization of clinical trials for promising LAMA2 CMD treatments.

The clinical and genetic spectrum of autosomal-recessive TOR1A-related disorders.

In the field of rare diseases, progress in molecular diagnostics led to the recognition that variants linked to autosomal-dominant neurodegenerative diseases of later onset can, in the context of biallelic inheritance, cause devastating neurodevelopmental disorders and infantile or childhood-onset neurodegeneration. TOR1A-associated arthrogryposis multiplex congenita 5 (AMC5) is a rare neurodevelopmental disorder arising from biallelic variants in TOR1A, a gene that in the heterozygous state is associated with torsion dystonia-1 (DYT1 or DYT-TOR1A), an early-onset dystonia with reduced penetrance. While 15 individuals with AMC5-TOR1A have been reported (less than 10 in detail), a systematic investigation of the full disease-associated spectrum has not been conducted. Here, we assess the clinical, radiological and molecular characteristics of 57 individuals from 40 families with biallelic variants in TOR1A. Median age at last follow-up was 3 years (0-24 years). Most individuals presented with severe congenital flexion contractures (95%) and variable developmental delay (79%). Motor symptoms were reported in 79% and included lower limb spasticity and pyramidal signs, as well as gait disturbances. Facial dysmorphism was an integral part of the phenotype, with key features being a broad/full nasal tip, narrowing of the forehead and full cheeks. Analysis of disease-associated manifestations delineated a phenotypic spectrum ranging from normal cognition and mild gait disturbance to congenital arthrogryposis, global developmental delay, intellectual disability, absent speech and inability to walk. In a subset, the presentation was consistent with foetal akinesia deformation sequence with severe intrauterine abnormalities. Survival was 71%, with higher mortality in males. Death occurred at a median age of 1.2 months (1 week-9 years), due to respiratory failure, cardiac arrest or sepsis. Analysis of brain MRI studies identified non-specific neuroimaging features, including a hypoplastic corpus callosum (72%), foci of signal abnormality in the subcortical and periventricular white matter (55%), diffuse white matter volume loss (45%), mega cisterna magna (36%) and arachnoid cysts (27%). The molecular spectrum included 22 distinct variants, defining a mutational hotspot in the C-terminal domain of the Torsin-1A protein. Genotype-phenotype analysis revealed an association of missense variants in the 3-helix bundle domain to an attenuated phenotype, while missense variants near the Walker A/B motif as well as biallelic truncating variants were linked to early death. In summary, this systematic cross-sectional analysis of a large cohort of individuals with biallelic TOR1A variants across a wide age-range delineates the clinical and genetic spectrum of TOR1A-related autosomal-recessive disease and highlights potential predictors for disease severity and survival.

Flavones provide resistance to DUX4-induced toxicity via an mTor-independent mechanism.

Facioscapulohumeral muscular dystrophy (FSHD) is among the most common of the muscular dystrophies, affecting nearly 1 in 8000 individuals, and is a cause of profound disability. Genetically, FSHD is linked to the contraction and/or epigenetic de-repression of the D4Z4 repeat array on chromosome 4, thereby allowing expression of the DUX4 gene in skeletal muscle. If the DUX4 transcript incorporates a stabilizing polyadenylation site the myotoxic DUX4 protein will be synthesized, resulting in muscle wasting. The mechanism of toxicity remains unclear, as many DUX4-induced cytopathologies have been described, however cell death does primarily occur through caspase 3/7-dependent apoptosis. To date, most FSHD therapeutic development has focused on molecular methods targeting DUX4 expression or the DUX4 transcript, while therapies targeting processes downstream of DUX4 activity have received less attention. Several studies have demonstrated that inhibition of multiple signal transduction pathways can ameliorate DUX4-induced toxicity, and thus compounds targeting these pathways have the potential to be developed into FSHD therapeutics. To this end, we have screened a group of small molecules curated based on their reported activity in relevant pathways and/or structural relationships with known toxicity-modulating molecules. We have identified a panel of five compounds that function downstream of DUX4 activity to inhibit DUX4-induced toxicity. Unexpectedly, this effect was mediated through an mTor-independent mechanism that preserved expression of ULK1 and correlated with an increase in a marker of active cellular autophagy. This identifies these flavones as compounds of interest for therapeutic development, and potentially identifies the autophagy pathway as a target for therapeutics.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Neuromuscular disorders: finding the missing genetic diagnoses.

Neuromuscular disorders (NMDs) are a wide-ranging group of diseases that seriously affect the quality of life of affected individuals. The development of next-generation sequencing revolutionized the diagnosis of NMD, enabling the discovery of hundreds of NMD genes and many more pathogenic variants. However, the diagnostic yield of genetic testing in NMD cohorts remains incomplete, indicating a large number of genetic diagnoses are not identified through current methods. Fortunately, recent advancements in sequencing technologies, analytical tools, and high-throughput functional screening provide an opportunity to circumvent current challenges. Here, we discuss reasons for missing genetic diagnoses in NMD, how emerging technologies and tools can overcome these hurdles, and examine future approaches to improving diagnostic yields in NMD.

Decoding the genetics of rare disease: an interview with Monkol Lek.

Monkol Lek, Assistant Professor at Yale University School of Medicine, and Associate Editor at Disease Models & Mechanisms, dedicates his research to finding a genetic diagnosis and improving treatments for rare disease patients. As he originally studied computer engineering at the University of New South Wales in Sydney, Australia, he now utilises computational methods to optimise large-scale genetic studies, provide globally accessible resources for genetic research communities and, importantly, resolve diagnostic odysseys for rare disease patients. Monkol completed his PhD in Prof. Kathryn North's lab at the University of Sydney, studying the genetics of muscle strength and performance, and then continued his investigation of muscle disease in Prof. Daniel MacArthur's lab at Massachusetts General Hospital and the Broad Institute. During his postdoc, he led several large-scale studies aimed at distinguishing pathogenic from benign variants, including the Exome Aggregation Consortium (ExAC) project ( Lek et al., 2016). Monkol established his own lab at Yale University School of Medicine, which continues to improve the diagnosis and treatment of rare muscle disease, and also focuses on underserved populations, whose genetic mutations are not as well characterised as those of European ancestry. In this interview, Monkol discusses how his own diagnosis with limb girdle muscular dystrophy has shaped his career and what he envisions for the future of genetic research in rare disease.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function.

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

MitoVisualize: a resource for analysis of variants in human mitochondrial RNAs and DNA.

SUMMARY: We present MitoVisualize, a new tool for analysis of the human mitochondrial DNA (mtDNA). MitoVisualize enables visualization of: (i) the position and effect of variants in mitochondrial transfer RNA and ribosomal RNA secondary structures alongside curated variant annotations, (ii) data across RNA structures, such as to show all positions with disease-associated variants or with post-transcriptional modifications and (iii) the position of a base, gene or region in the circular mtDNA map, such as to show the location of a large deletion. All visualizations can be easily downloaded as figures for reuse. MitoVisualize can be useful for anyone interested in exploring mtDNA variation, though is designed to facilitate mtDNA variant interpretation in particular. AVAILABILITY AND IMPLEMENTATION: MitoVisualize can be accessed via https://www.mitovisualize.org/. The source code is available at https://github.com/leklab/mito_visualize/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.